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Abstract We investigate some peculiarities in the calculation of the two-loop beta- 
function of A = 1 supersymmetric models which are intimately related to the so- 
called “Anomaly Puzzle”. There is an apparent paradox when the computation is 
performed in the framework of the covariant derivative background field method. 
In this formalism, it is obtained a finite two-loop effective action, although a non¬ 
null coefficient for the beta-function is achieved by means of the renormalized 
two-point function in the background field. We show that if the standard back¬ 
ground field method is used, this two-point function has a divergent part which 
allows for the calculation of the beta-function via the renormalization constants, 
as usual. Therefore, we conjecture that this paradox has its origin in the covariant 
supergraph formalism itself, possibly being an artifact of the rescaling anomaly. 

Keywords Renormalization Regularization and Renormalons • Superspaces • 
Renormalization Group • Gauge Symmetry 


1 Introduction 

In particle physics, symmetries have always been used as guides in order to con¬ 
struct theories to describe nature, an idea that culminated in the Standard Model 
(SM) itself. Although it has passed many experimental tests, the SM must be 
viewed as an effective theory since it does not incorporate all fundamental in¬ 
teractions. Therefore, extensions to it have been proposed, supersymmetry being 
one of the most appealing from the theoretical viewpoint. The reason lies on the 
elegance of its construction since it is a natural extension of the Poincare group. 
Much effort have been dedicated to the subject after it was first proposed in the 
seventies HEE]. 
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The investigation of quantum corrections to supersymmetric models has very 
interesting peculiarities. For example, the choice of an adequate regularization 
technique is a highly nontrivial subject. A model which incorporates supersym¬ 
metry is dimension specific, the spacetime dimension playing an important role in 
the matching between the bosonic and fermionic degrees of freedom. This poses 
restrictions in the use of methods based on the analytical continuation of the 
spacetime dimension. Regardless the method of calculation, the quantization of 
supersymmetric models presented some delicate questions, such as the so called 
“Anomaly Puzzle”, which can be summarized as follows. 

Since the works of Piguet and Ferrara [UOinilT] , it is known that there exist two 
supermultiplets, one for classically conserved supercurrents and the other a chiral 
supermultiplet for the scale anomalies. In the first multiplet are the classically 
conserved currents associated to the U(l) chiral R invariance, to supersymmetry 
and to translation invariance, with the last two conserved at all orders, while the 
R current is not. In the chiral supermultiplet are the scale anomalies associated 
with the R current, the trace of the supersymmetry current and the trace of the 
energy-momentum tensor. 

At the core of the Anomaly Puzzle is that an unique axial current operator only 
exists at tree level. At quantum level the R current in the supermultiplet is broken 
with a coefficient which is proportional to the beta function of the scale anomaly, 
which can have corrections to all orders. However, according to the Adler-Bardeen 
theorem [8], the chiral anomaly is exhausted at one-loop order in perturbation 
theory. Remarkably, its coefficient is given by the one-loop value of the gauge beta 
function. Despite these differences in the quantum realizations of the axial current, 
Piguet and collaborators obtained a relation which links the beta function of the 
scale anomaly to the nonrenormalized coefficient of the axial current anomaly 
ElllQ]. Since the work of Novikov, Shifman, Vainshtein e Zakharov (NSVZ) |11| . 
which obtained an exact expression for the beta function of V = 1 Super Yang- 
Mills (SYM) theory, many other works |12lll5l[T8lll91l21l[^l23ll241l251l26| followed 
in which different regularization methods were applied and in all cases higher order 
corrections for the beta functions were found. 

This controversial result attracted great attention and different explanations 
were provided. According to Shifman and Vainshtein m, it is necessary to dis¬ 
tinguish between the Wilson effective action and the sum of vacuum loops in the 
external fields. The first renormalizes only at one-loop level whereas the second 
receives higher order contributions due to infrared modes. On the other hand, 
Arkani-Hamed and Murayama m argued that the solution to the problem can be 
stated in a way independent of the infrared modes using the distinction between 
the holomorphic gauge coupling and the canonical gauge coupling. According to 
the authors, the dilatation anomaly is in the same multiplet of the anomaly 

and is exact at one-loop order. However, due to the anomaly, the vectorial mul¬ 
tiplet does not possess canonical kinetic terms after the dilatation. In order to 
get canonical kinetic terms in the vectorial multiplet an additional change in the 
normalization is needed. Therefore, the anomaly coming from the modified dilata¬ 
tion is not in the same multiplet of the anomaly and receives contributions 

beyond one-loop order. This argument was criticized in some papers, since to keep 
the low energy physics unchanged, it is necessary to take into account the infrared 
modes in the derivation of the anomaly. In this sense, it is somewhat equivalent 
to consider the scale anomaly or to calculate the expectation value of the Wilson 
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effective action m- In [16], it was claimed that since the definition of the gauge 
coupling of A'" = 1 SYM may depend on the renormalization scheme, so does the 
beta function. They showed that the trace anomaly is one loop exact in a certain 
scheme, the important point being to examine in which scheme the quantum action 
principle is valid. To summarize, as observed in m, although the i?-current and 
the stress tensor belong to the same classical supercurrent, in the quantum regime 
it bifurcates. It is not possible to construct an unique quantum supermultiplet 
which contains both the stress tensor and the i?-current. 

This discussion above also appears in a perturbative analysis. Within the su¬ 
pergraph approach to supersymmetric models, along with on-shell infrared di¬ 
vergences of Yang-Mills theory, additional off-shell infrared divergences appear 
which must be distinguished from ultraviolet ones before renormalization is car¬ 
ried out. The mixing of these two types of divergences is in the center of this 
debate. A consistent approach should proportionate an unambiguous distinction 
between the infinities involved and the arbitrary scales which are byproducts of 
the subtractions. In dimension-type regularizations mM the two-loop correction 
to the /3-function comes from a local evanescent operator, which would be absent 
in the physical spacetime dimension. So, Grisaru, Milewski and Zanon conjectured 
that no divergence should occur beyond one loop. This is true, as we will see, de¬ 
pending on the approach adopted in the calculation. However, even in the case 
where the divergences do not occur beyond one-loop order, this does not mean 
the two-loop /3-function vanishes. Instead, the derivation of the renormalization 
group functions needs some reinterpretation, which appears to be related to scal¬ 
ing anomaly [l3]- In four spacetime dimensions. Differential Renormalization was 
applied m in the evaluation of the two-loop /3 function. It was found that the 
result depended on infrared modes, which play a passive role. Moreover, within 
differential renormalization, it was found that the scale referring to one-loop renor¬ 
malization is the one to give rise to the two-loop coefficient. Because differential 
renormalization delivers finite renormalized amplitudes by construction, it would 
be interesting to investigate how renormalization is effected within an invariant 
framework which both operates in the physical dimension and displays explicitly 
the ultraviolet behavior in terms of the renormalization constants. In |26j, a four 
dimensional regularization framework was used in the computation of the two- 
loop coefficient of the SYM beta-function with the use of the background field 
method in the covariant derivative formalism [28] . Due to the non-abelian charac¬ 
ter of SYM, the background field method is urged to be applied, since it results 
in a huge simplification in the number of diagrams. It was found that there is no 
two-loop divergence which, in a first view, could indicate the absence of higher 
loop corrections to the beta-function. However, the renormalized two-point Green 
function still depends on the renormalization scale introduced at one-loop level, 
allowing the computation of the two-loop coefficient for the beta function, which 
was shown to be non-null. It is interesting to explore such result also in view of the 
property that, if the n-|-1 loop coefficient of the beta function of a supersymmetric 
theory vanishes, then it is finite to n loops 

Finally, it is interesting to discuss the dependence of the /3 functions in terms 
of the renormalization scheme applied [16] . As discussed in El, the holomorphic 
beta function is one loop exact whereas the canonical is given by the NSVZ /3 
function. The canonical gauge coupling comes from the canonical normalization of 
the holomorphic gauge superfields which is anomalous and is determined by the 
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axial anomaly. Therefore, while the holomorphic coupling defines the Lagrangian 
(for notation see please next section) 

Ch = ^ j d^ew%Vh)w%Vh) + h.c. ( 1 ) 

the canonical coupling defines Cc by replacing l/(/^ with 1/g^ — idyM /as well 
as V/i —>■ gVc which however are not equivalent. Because of such a rescaling anomaly 
the (holomorphic and canonical) coupling constants turn out to be related by a 
non-local relation. 

In this context, several renormalization scheme dependence issues may arise. 
For example, for A'" = 1 massless supersymmetric QED (SQED) regularized by 
High Covariant Derivatives the NSVZ /3-function is naturally obtained for 
renormalization group functions defined in terms of the bare coupling constant and 
do not depend on the renormalization prescription. However, if defined in terms 
of the renormalized constant, the NSVZ /3-function is only obtained in a special 
subtraction scheme namely the NSVZ scheme. Thus, for the exact /3-function it is 
natural to ask in which (precise) scheme its expression holds. For instance the /3- 
function in the minimal subtraction scheme of dimensional reduction is not given 
by the NSVZ /3 function beyond two loop level [29]. In any case, it should be 
noticed that renormalization scheme dependence arises when comparing different 
approaches. 

In this work, motivated by our results in [26], we use massless SQED as a 
laboratory to investigate peculiarities in the calculation of beta-functions of su¬ 
persymmetric theories up to two-loop order. Particularly, we would like to under¬ 
stand the apparent paradox found in |26| in the simplest context. Therefore, in the 
present work we will use SQED as a probe. For this purpose, we will compute the 
two-loop coefficient of the SQED beta-function using two different approaches: 
the standard background field method m and the one based on the covariant 
derivative formalism [28]. We will find that in the first case there is a two-loop 
divergence, allowing the computation of the beta-function coefficient by standard 
renormalization constants. It is also possible to perform the computation using the 
renormalized two-point Green function, furnishing the same result as before, as 
expected. In the second case, we obtain the same result of [26]: there is no two-loop 
divergence, even though the renormalized two-point Green function depends on a 
renormalization scale, furnishing a non-null value for the two-loop beta-function 
coefficient. Therefore, we conjecture that the paradox found in [2^ has its ori¬ 
gin in the covariant supergraph formalism itself, possibly being an artifact of the 
rescaling anomaly, discussed in [T1[T1[2Z1[13]. Indeed, a mechanism of corrections 
to the one-loop result from one-loop anomalies is described in El, through the 
quantum breaking of holomorphy of the coupling constant. 

The outline of this paper is as follows. In section H, we present the super- 
symmetric QED and evaluate, in the framework of the standard background field 
method, the two-loop /3-function, using both the renormalization constants and 
the renormalized two-point function in the background field. In section HI, the 
formalism of covariant derivative background field method is applied in the cal¬ 
culation of the two-loop beta-function of SQED. We present in section IV our 
conclusions and perspectives and some results of integrals are displayed in the 
appendix of section V. 
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2 N = 1 SQED in the standard background field method 

In the superfield formalism, the classical action of the massless = 1 supersym¬ 
metric quantum electrodynamics (SQED) is given by [31] 

S = I d‘^xcPew‘^+ j d!^xd!^e§+e^^$+ + I d'^xct^e (2) 

where ^ is a chiral field that express the matter part of the action and E is a real 
scalar superfield that contains the gauge field A^i of QED as one of its components 
(therefore, it is the supersymmetric generalization to the gauge field). Finally, W 
is the supersymmetric generalization to the stress tensor of QED. In terms of the 
superfield V, one has 

J d'^xd^e W'^ = ^ j d‘^xd‘^6 VD'^D'^DpV. (3) 

The following step would be to perform the quantization of the classical theory. 
However, since we want to use the background field method, we have to introduce 
this new field at this point. For the abelian case, we will have a linear quantum- 
background splitting as below |32| 


V^V + B, (4) 

where B is the background gauge field. We may now perform the quantization as 
usual, introducing a gauge-fixing term for the quantum gauge field E, as well as 
sources [30]. The relevant fact to be noticed is that, by construction, the action 
will be gauge invariant in the background gauge field. This must remain valid even 
after renormalization. Thus, the renormalization constant for the background field 
Zb will be related to the one for the gauge coupling Zg as follow^ 

ZgZl^^ = 1. (5) 

Thus, in order to obtain the beta-function of the theory, we need only to compute 
the two-point functions with background fields as external legs. 

As stated in the introduction, we intend to compute the two-loop corrections 
of the SQED beta-function. The Feynman rules can be derived from the action 
[32| , and the relevant ones for our computation are expressed in figure [T| 

We start by the one-loop contribution, whose diagrairQ, using the background 
field method, is depicted in figure (H furnishing the following effective action 

yl(l)=2^/’ Bi-p,e) [ 

^ Jp,e Jk 

^We are using the definitions go = Zgg and Bq = where ^'q, Bq and g, B are bare 

and renormalized functions, respectively. 

^We do not include a tadpole diagram since, as we are dealing with a massless theory, it can be 
promptly set to zero in the Implicit Regularization formalism. However, even if such diagram 
was included, it would cancel out in the sum. Thus, in order to simplify the discussion, we opt 
to omit hereafter all the tadpole diagrams. 


- i 


DjjDg - fc 


k?{k p)2 


B{p,e), 


( 6 ) 
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Fig. 1 Feynman rules needed for the evaluation of two point functions in the background field 
up to two-loop order. 


where we have already performed the D-algebra manipulations, stands for 
f g for f f 9 and B{p,6) is the background gauge field. The fac¬ 

tor 2 accounts contributions from chiral helds <P with different signs. Regarding 
supersymmetric definitions and conventions, we are following the ones found in 

El. 



Fig. 2 1-loop diagram. 


To proceed, we need to resort to some regularization technique. We will choose 
the Implicit Regularization (IReg) formalism |331I34| . which, by not resorting to 
any kind of dimensional extension can be promptly applied in supersymmetric 
theories |351I361[^ . The method resorts recursively to the mathematical identity, 

_ 1 ^ 1 _ Pi +‘i‘Pi-k _ 

{pi -\- k)'^ — m? k? — (fe^ — m^) [(pj -j- — m^] ’ 


( 7 ) 
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in order to extract the external momenta, pi, from the divergent integrals. In the 
case of massless models, a fictitious mass, , is used, which is eliminated from 
the result by means of the limit ^ 0 and of scale relations. As a byproduct, a 
mass scale is introduced, which is adequate for the computation of the renor¬ 
malization group functions. A detailed account on IReg can be found in |361I371 
[38]. 

After some D-algebra manipulation the final result is 


B{-p,9)D^D^DpB{p,9) 


i+i)^ J^^B{-p,9) + 2fe2j r^^’^^B{p,9), 


Ii^g{X^)-b\n[-^]+2b 


+ 


( 8 ) 


with b = i/(47r)^, 


4s / (fc2 - /r2)2 ' ( A2^ ) 

and 

gtJ-l'- tJ-j = / 

Jk 

where ■a'j = ... gk-j-ik-j -|-symmetric combinations and the index i in¬ 
dicates the superhcial degree of divergence of the surface term. In we 

omitted the upper index for Z = 1. 

Some comments are in order: notice that, apart from the surface term Rq ’ , 
we have a gauge invariant result. This was expected, since, as we are working 
with the background held method, gauge invariance is explicitly maintained being 
broken just by regularization dependent terms. Thus we verify that the condition to 
preserve gauge invariance, even in supersymmetric theories, is to set surface terms 
to zero, as discussed in (SB] . From now on, we will not display the one-loop or higher 
order surface terms, which will all be set to zero. Another aspect to be mentioned is 
the appearance of a divergent integral parametrized as a //og(A^). Such term must 
be renormalized as usual and it will contribute to the renormalization constant 
Zb, as we are going to show in the end of this section. Notice also the appearance 
of the renormalization scale A^ in the hnite part as well. 

We proceed now to the two-loop contributions, which are depicted in hgurejSj 
The hrst diagram furnishes, 


d 


fc'"" • • • 


dfc, 


hi (A:2 -p2)- 


■Ih 


l-i 


A2 


( 10 ) 






/ B{-p,9) / 

Jp,9 Jk 

fc2(fc — p)2 

(ifo) 



+ 


2^ / B{-p,9) I 

^ Jp,e Jk 


-h 

k?(k — p)2 


B{p,9) 


{-i)Iiog[X^)B{p,9), (11) 


where the second line is a subdivergence which is be to canceled out by a one-loop 
counterterm according to the subtractions stated by the Bogoliubov’s recursion 



















formula. Instead of evaluating the first line, we proceed to the next diagram (the 
reason will be apparent shortly), which has the result quoted below, 


y^(2) _2i 
^102 - ^ 


/■ / ^ /■ 

- k^^DpDa - k^ 


\ { P^\ 1 

Jp,e Jk 

fc^(fc + p)2 

(ih) 



+ 2 - 


I [ 

J J k 


D'^D'^ -k^ 


k'^ik + p)^ 


B{p,e) 


{-i)llog{\^)B{p,e). ( 12 ) 


Once again, the last line will be subtracted by applying Bogoliubov’s recursion 
formula. 










Fig. 3 2-loop diagrams. 


The reason why we have not performed the computation of the integral in k will 
become apparent now. After performing the D-algebra on the diagrams 6i • ■ • 64 
we see that they all can be expressed in terms of A^ai and as below 


dp)=dp)=-A 


( 2 ) 




( 2 ) 


(13) 
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We proceed to diagram c, which is given by 

//‘-''’'''ipiyraryysfe'')- (») 

No further analysis should be taken, since the last diagram, d, can be written 
in terms of , and 

4"^ = 4"^ + 4"^ - 4"^ I B{-p,e)p^ [(/io 3 (A") + 6)4^’") - 

-2^ j^^B{-p,e) +2fer4’"))] B{p,e) 

-2^ j ^B{-p,e)D^D‘^DpB{p,e)b 7,„g(A2)-61n(^-^j -4^’") B{p,e). 

(15) 

Therefore, the final two-loop correction for the effective action, after discarding 
the surface terms, is given by 

A^^^ =-g^ j ^B{-p,e)D^D^Di^B{p,e)b Ii,g{X^)-b\n(^-^^ B{p,e). (16) 

A curious fact is that, although we are working at two-loop level, we have 
a divergent result parametrized by //og(A^), which is a typical basic divergent 
integral (BDI) of one-loop order. Our next task is to perform the renormalization 
of the theory, whose bare action is 

So = J d'^xd'^e ^BoD^D^DfiBo + <Po+e^°^<>'Po+ + 

-|- gauge fixing terms + terms on Vb- (17) 

Performing a multiplicative renormalization defined by Bq = Zg B, go = Zgg, 
and <Po± = Z^j_<P±, one can easily find that the counterterm for the two-point 
function in the background field is given by A = Zg — 1. As already mentioned, 
in the background field method, the relation ZgZ^ = 1 holds. Therefore, the 
beta function can be calculated through the two-point function renormalization 
constant as follows 

= = ( 18 ) 

Supposing that both the countertem A and the /3-function have expansions in the 
coupling constant g, we arrive at the expressions 

ft = (19) 

and 

ft = + 1^)1 - f ). (») 

where /3j is the i-loop coefficient of the beta function, Ai is the i-loop two-point 
function counterterm and A is the renormalization group scale. Using the one- 
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and two-loop corrections obtained in eqnations (|8]) and (HU), respectively, and 
adopting a minimal snbtraction scheme (which in the IReg framework amounts to 
the subtraction of BDI’s only), we have 

Ai = and Aa = 267j„g(A"). (21) 

Finally, by using 

X^^hogiX^) = 2X^^Ii,g{\^) = -2fe, 

with h = we obtain the contributions for the beta function 

two-loop level in the IReg formalism, 

^ ^ (4^®^ 8(47r2)2^^ + 

which agrees with previous ones found in the literature [ 39 ll 30 llin[nl FI 

It should be emphasized that we performed the above computation using only 
the renormalization constant for the background field, Zb- Alternatively, the com¬ 
putation can be performed using the renormalization group equation, in which 
only the renormalized effective action is considered. As a consistency check, we 
also compute the two-loop SQED beta-function using the renormalized effective 
action. 

As can be found, for example, in [12] , the renormalization group equation reads 


( 22 ) 

of SQED up to 


.9 

^dX '^^dg 


Giliig,X) = 0, 


(24) 


where 

9 = ^§^9 =-gX-^\nZg and j = X-^lnZ b- (25) 

Since in the background field method ZgZ^^^ = 1, we find that 7 = ^- Thus, 




(- 

\dg 



(g. A) = 0. 


(26) 


From equations (|8|) and m, we obtain the renormalized two-point function as 
below 


GiVn{g,X) = ^ j^^B{-p,e)D^D^DpB{p,e)i^l + ib 


-2 




(27) 


Replacing the expression above in equation (1241) . we finally obtain the same ex¬ 
pression as in equation (|23|). 

We finish this section with some comments. In [15], the authors conjectured 
that no divergence should occur beyond one-loop in A = 1 SYM theory if the cal¬ 
culation is performed in the physical dimension of the model. Here, we performed 
the two-loop calculation of the two-point function with the use of the standard 


®To obtain the results of the last three references, it should be taken into account that our 
definition for the coupling constant m differs from the usual one by a factor \/2- 
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background field method for SQED and found a divergent result. It should be ob¬ 
served, however, that the divergence found is typical of a one-loop calculation. In 
the next section, we will adopt an alternative approach, which is the background 
field method in the covariant supergraph formalism. It is just what was done in 
|2fi| for SYM. In this case we will see that the result will match with the conjecture 
of |45| , although this does not imply the two-loop beta-function is null. 


3 The SQED /3-function in the covariant supergraph formalism 

We now perform the whole calculation again, by using a different approach: the 
background field method based on the covariant supergraph formalism |28| . The 
reason is the following: in [26], in which the aforementioned formalism was applied, 
the Super Yang-Mills theory was studied and, particularly, the beta-function coef¬ 
ficients were computed up to two-loop order. It was found there that no two-loop 
divergence appeared,which could be taken as an indication of a null two-loop coef¬ 
ficient. However, the renormalized effective action at two-loop order still carried a 
dependence on the renormalization scale A, allowing the computation of the beta 
function from the renormalization group equation, furnishing a non-null result for 
the two-loop coefficient. Therefore, it seems that there is an inconsistency, since 
both approaches should be equivalent. It was conjectured there that this difference 
should have its origin in the rescaling anomaly, as suggested by [l3] in which the 
author discusses that in a framework that uses the canonical coupling (as ours), a 
modification of the usual multiplicative renormalization program should be neces¬ 
sary. Therefore we will in the following use SQED as a probe to study if the same 
behavior occurs in this case. 

A complete description of the background field method based on the covariant 
supergraph formalism can be found in |281I32| . The main idea is to take a step 
backward and work, from the beginning, with an action that depends only on 
background covariant derivatives. This way, all dependence on the background 
field will only appear implicitly. The main gain on this approach is the reduction 
in the number of diagrams. For instance, the two-loop correction we are going 
to compute requires considering only three diagrams, instead of eight, as in the 
previous formalism. Explicitly, the quadratic part of the action in the gauge fields 
we are going to work with is given by 



(28) 


where V is a covariant derivative in the unsplit gauge field {V + B). The splitting 
can be carried out, in the quantum-chiral but background-vector representation, 
as 



(29) 


being V background covariant derivatives. The quantization procedure should be 
carried out from this point, adding gauge fixing and source terms, as usual. Chiral 
fields should also be included to define SQED properly (this fields must also be 
written in the background covariant representation). After all these considerations. 
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we obtain the covariant Feynman rules [32] which, applied to our case, furnish the 
following one-loop effective action (the diagram depicted is the same of fig. |2|) 

= (30) 


Notice that our result has no explicit dependence on the background held B. It 
appears only through the held strength W“ and the spinor connection Fa- This is 
a feature of the method, since background covariance, by construction, is always 
maintained. We obtain: 


= (- 


(-*5^) / W“(p)ra(-p) 

Jp,e 


7,„g(A")-61n(-^) +2fe 


(31) 


To conclude the one-loop calculation, we write our result in terms of an explicit 
background gauge held. For this purpose, we recur to the dehnitions of W“ and 
Fa found in |31] . 

Fa = S-Ild Wa = iD^DaB, (32) 

such that our hnal one-loop result reads 





B{-p,e)D^D^DpBip,e) 





(33) 


Comparing with our previous expression, eq. (|8|), we notice that there is no 
dependence on a surface term this time. We conjecture that this feature may be 
a consequence of the method which, by maintaining background covariance from 
the beginning, have automatically canceled all gauge-breaking terms that could 
occur. We also remark that this is the only difference between the two results. 

We proceed now to the two-loop contribution whose diagrams are depicted in 
hgure m 





(a) 


(b) 


(c) 


Fig. 4 2-loop diagrams in the background covariant approach. 


As in the one-loop case, the effective action will have only an implicit depen¬ 
dence on the background gauge field, through the field strength W“ and the vector 
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connection F—, as the expressions for the diagrams above reveal, 




1 . 


4^^ = 4/ 


Jp,e L 


/ [ 
Jp,9 . 


2W“(-p,e)W“(p,e) (7i)^^ + -V“Wa(-p,0)vPw^(p,e)72 


--r^-p,e)rHp,e) (h)^ 


-r^-p,9)ra(p,e)h 


(34) 

(35) 


and 


= 4/ / \jr^{-p,e)rHp,e) {h)^, + iv“w„(-p,0)v^w^(p,0)76 
Jp,e ~ ^ 


where 7^ are the following integrals 


, k + k)^k^(k — pY ’ 


(7i) .=(7^.(7 i) =a>^.{-Af 

V ^ J cxa act v aaV / / 

«/OJ 

h = (-*)" / 

«/ aJ 


, fc <J^('3' + kYk‘^{k — pY 


(36) 

(37) 

(38) 


(■^ 3 )^ = (^ 3 )^,, = [ 

Jq 

I,k qHp + kYk^' 

i-if f 

J a. 


Akfj,kii ‘2p^kif ‘2k^pif p^pu /on\ 

,,k qY9 + kYk\k-pY ’ ^ ^ 


74 = — 


ih)a, = 


_ p V 


ab ^aa^^/3 ^aa^^/3 


_ fl V 


^ ^.A.^ a qK 


4fc^gt/ - 2k^iPi, + 2pfj,qu - p^pv 

k q‘^{q + kYk‘^{k + pY{q-pY 


h = (- 


{-A / 

«/ a. 


k q‘^iq + kYk‘^{k+pY{q-pY' 


(40) 

, ( 41 ) 


(42) 


Since now we have some integrals with off-shell infrared divergences, we will 
explain with same detail the treatment of the first integral. Within the IReg pro¬ 
cedures, we have 




(p-fc)l 


k {k^-pYAk-PY-A] [ 


(X^)-bin I- 


[k^-pY 

A2 


+ 2fo 


(43) 


where we have to notice the inclusion of the fictitious mass which must be 
added in order to regularize the infrared divergence. There is also an UV divergence 
parametrized as a Iiog{A) which is just an one-loop subdivergence that is going 
to be canceled by the application of Bogoliubov’s recursion formula. After the 
subtraction of the subdivergence one obtains 


= i-A [pm-F + - 2hU^ 


with 


f/M = 


Ik k^{k-pY' 


uA = 


Jk k^{k-p)‘- 


• In — 


A2 


(44) 

(45) 

(46) 
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and 




L ’ _ 1 1 

“ ^^0./, (fe2-p2)2p_p)2_^2] 


—6 In 


A2 j 


+ 2b 


(47) 


As we can see, the infrared divergence is concentrated in the F integral. To proper 
treat the IR divergence, one could, for example, resort to the IReg generalization 
presented in |46| . However, in the present case the integral F will cancel ont with 
other contributions, not reqniring any further treatment. For the other integrals, 
we obtain 


h (48) 

={-if [pt^P^^F + 2b{p^,ui^'> - 2pf,Uu+p^U^^^ - 2p^Ut, - 2uj,V + 4(7^,.)] , 

(49) 


h ={-if 


log' 

— 1/1 7^2 


-bll2i^^)+2bIi,g{X^) + ^-\n^(-^] -b^\n(-^] +p^F 


A2 




h=r 


o 


(50) 

(51) 

(52) 


where we defined 


Unu = 


L 


kuky 


k k^{k-pY' 


1°^ 


I^ = 
^11 — 


I^ = 

J-IJ — 


r02 = 
^ fll/ - 


f kf,K , ( k^\ 

J^k^{k-py'^\ 

I 
I 
I 
I 


,fc ( 72 (g + fc)2fc2(fc+p)2(g-p)2’ 
kfj, 

,/c q^{q + k)^k^{k + py{q-py' 
_ <F _ 

,/c q^{q + k)'^k'^{k + py{q-p)'^’ 
kfi^y 

,/c q^{q + kYk'^{k + pY{q-py 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


The results of the integrals can be found in the appendix. We now proceed 
noticing that, in the effective actions we have different structures in terms 

of the field strength and vector connection. For reasons that are going to be ap¬ 
parent soon, we choose to group all the contributions proportional to the vector 
connection, obtaining 




4/ ^ ^ r^{-P, e)rHp, o) {h)^ + i {g^h) + i {h)^ 


(59) 
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Replacing the values of the integrals found in the appendix, we have, after 
discarding the surface terms 




f r-i 

J 


i-p,e)rHp,e) 


PaPb 


n2 3^ 


6 ^1 n(-^ 


+ 


0 TV 


36 


+ 


b^C(3) ^ Fp 


(60) 


Notice that, although we have dealt with ultraviolet and infrared divergent 
integrals, the net result is hnite and gauge invariant, and obeys the following 
relation 


Id^er^{-p,e)r^{p,e) = ^Jd^ew^Wc.. (6i) 

Since we also have the relations, 

I dW“(-p,0)p„<iW“(p,e) = iy d^ew‘^{-p,e)Wc.{p,e) (62) 

and 


I d^dv“Wa(-p,d)v^w^(p,d) = -^ Id‘^ew°‘{-p,e)Waip,e), (63) 

and (A)c,q, ocpaa, we finally obtain the two-loop effective action. 


^(2) ^4^4 f d2gW“(-p,0)Wa(p,d)6^ 

Jv 


^In 


A2 


3C(3) 


24 


-2 


(64) 


which, written in terms of the background held by means of eq. (El, is given by 


yl( 2 ) = ^ f B{-p,9)D^D^DpB{p, 
^ Jp,6 


e)b^ 


21n(-f2 ) +- + 6C(3)-8 


■ (65) 


Some comments are in order. One should compare the result above with the 
one expressed by eq. (I16II which was obtained using the standard background held 
method. A notorious difference is the disappearance of the UV divergent inte¬ 
gral. This means that there isn’t a two-loop contribution for the renormalization 
constant Zb, which, at a hrst view, could lead one to think that the two-loop 
coefficient of the SQED beta-function is null. However, one could also compute 
the beta-function using the renormalization group equation. For this purpose, the 
following renormalized two-point function is needed 




l + ib 


In - 


A2 


-2 


+2b-‘ 


In - 


A2 


+ - + 3C(3)-4 




( 66 ) 


By comparison with eq. EH). one notices that in both methods (standard and 
covariant derivative background held method) the dependence of the renormalized 
two-point function in the renormalization scale A is the same. Since this is the only 
relevant part for the computation of the SQED beta-function, we hnd that in both 
cases there is a non-null two-loop beta-function coefficient given by 


P2 


8(47r2)2 


(67) 
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4 Discussion of the results and perspectives 

In this paper we have studied massless SQED up to two-loop order. Our purpose 
was to study the intriguing fact that different approaches in the use of the back¬ 
ground field method result in the existence or not of a divergent part in two-loop 
calculations in A'" = 1 supersymmetric theories, even though the corresponding 
/3-function coefficient is the same. We used the Implicit Regularization framework, 
since it operates in the physical dimension of the theory (respecting supersymme¬ 
try) as well as displays in a clear way UV and IR divergences and regularization 
dependent surface terms. The use of the background field method simplihes con¬ 
siderably the calculations by reducing the computation of the beta-function to the 
knowledge of two-point functions in the background field. The two approaches used 
were the following: the standard |30] and the covariant derivative [28] background 
held method. In the Hrst case, we obtained that the one and two-loop effective 
action contained a divergence. Therefore, the beta-function could be computed in 
the usual way, by dehning a renormalization constant in the background held. On 
the other hand, by using the covariant derivative background held method, we 
obtained that the two-loop effective action had no divergence. This could imply 
that the beta-function would not receive higher order corrections. However, the 
renormalized two-point function still depended on the mass scale A, which allowed 
us to obtain the two-loop /3-function coefficient. Both approaches yielded the same 
result 

coinciding with the one obtained before in the literature |llll391linil41| . It should 
be noticed that even in the case of the standard background method the beta- 
function could be computed by using the renormalization group equation, which 
delivered the same result as before. 

We emphasize that our main point is to find out if the computation via the 
standard or covariant derivative background held method could, respectively, give 
rise or not to the explicit divergent behavior at two loop order, yet there is no 
doubt about the value of the two loop correction to the beta function as they 
agree (computing via the renormalization constants or via RG equation). In other 
words, there is no doubt about renormalization scheme in our analysis, since the 
value of the beta function obtained in both methods coincide with each other 
corroborating the universality of the two loop coefficients of the /3-function. The 
question would be in which method the multiplicative renormalization program is 
still applicable, since the divergent behavior of the effective action in both methods 
is not the same. 

This particularity for the calculation with the covariant derivative background 
held method was already obtained in the context of SYM theory [26]. Therefore, 
we found out with our computation that the above behavior is not characteristic 
of the SYM theory, being shared by the SQED theory as well. This allows us 
to conjecture that the reason may lie on the rescaling anomaly and the usual 
multiplicative renormalization program should be modihed as suggested by |43| . 
being inherent to the dehnition of the covariant derivative background held method 
itself. 

As perspectives we should include the study of how exactly the rescaling 
anomaly manifests itself in the covariant derivative background held method. Thus, 
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one expects to be able to introduce some modifications in the usual multiplicative 
renormalization in order to solve this controversy in the computation of the beta 
function. 


Appendix A: List of integrals used in this work 

For the integrals needed in section [3l we have the following results: 
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— 


Uuu = 


^ av — 
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(A.5) 

(A.6) 

(A.7) 




— —Oki^ { -^ ( “ 


(A.8) 


Acknowledgements A. Cherchiglia acknowledges fruitful discussions with M. Perez-Victoria 
and thanks Universidad de Granada for the kind hospitality. A. Cherchiglia and M. Sampaio 
acknowledge financial support by FAPEMIG and CNPq, Conselho Nacional de Desenvolvi- 
mento Cienti'fico e Tecnologico - Brazil. 


References 

1. J.-L. Gervais and B. Sakita, Nucl. Phys. B34, 632 (1971). 

2. D. Volkov and V. Akulov, Phys. Lett. B46, 109 (1973). 

3. P. Ramond, Phys. Rev. D3, 2415 (1971). 

4. S. Ferrara and B. Zumino, Nucl. Phys. B87, 207 (1975). 

5. T. Clark, O. Piguet, and K. Sibold, Nucl.Phys.B143, 445 (1978). 

6. O. Piguet and K. Sibold, Nucl. Phys. B196, 428 (1982). 

7. O. Piguet and K. Sibold, Nucl. Phys. B196, 447 (1982). 

8. S. L. Adler and W. A. Bardeen, Phys. Rev. 182, 1517 (1969). 

9. O. Piguet and K. Sibold, Int. J. Mod. Phys. A1 (1986) 913. 

10. C. Lucchesi, O. Piguet and K. Sibold, Helv. Phys. Acta 61 (1988) 321. 















18 


11. V. Novikov, M. A. Shifman, A. Vainshtein, and V. I. Zakharov, Phys. Lett. B166, 329 
(1986). 

12. W. Siegel, Phys. Lett. B84, 193 (1979). 

13. M. A. Shifman, A. I. Vainshtein, Nucl. Phys. B277, 456 (1986). 

14. Nima Arkani-Hamed, Hitoshi Murayama, JHEP 0006, 030 (2000). 

15. J. Mas, M. Perez-Victoria, and C. Seijas, JHEP 0203, 049 (2002), arXiv:hep-th/0202082 
[hep-th]. 

16. K. Yonekura, JHEP 1203, 029 (2012). 

17. Xing Huang, Leonard Parker, Eur. Phys. J. C71, 1570 (2011). 

18. K. V. Stepanyantz, Theor. Math. Phys. 142, 29 (2005). 

19. K. V. Stepanyantz, Nucl. Phys. B852, 71 (2011). 

20. K. V. Stepanyantz, JHEP 1408 (2014) 096. 

21. L. Avdeev, O. Tarasov, and A. Vladimirov, Phys. Lett. B96, 94 (1980). 

22. M. T. Grisaru, M. Rocek, and W. Siegel, Phys. Rev. Lett. 45, 1063 (1980). 

23. W. E. Caswell and D. Zanon, Phys. Lett. BlOO, 152 (1981). 

24. A.Pimenov,E.Shevtsova, andK.Stepanyantz, Phys. Lett. B686, 293 (2010), 

arXiv:0912.5191 [hep-th]. 

25. E. Abdalla and R. Jasinschi, Nucl. Phys. B286, 42 (1987). 

26. H. Fargnoli, B. Hiller, A. B. Scarpelli, M. Sampaio, and M. Nemes, Eur. Phys. J. C71, 
1633 (2011), arXiv:1009.2976 [hep-th]. 

27. E. Kraus, Nucl. Phys. B620, 55 (2002). 

28. M. T. Grisaru and D. Zanon, Nucl. Phys. B252, 578 (1985). 

29. 1. Jack, D. R. T. Jones and C. G. North, Nucl.Phys.B 486 (1997) 479. 

30. L. Abbott, Nucl. Phys. B185, 189 (1981). 

31. J. Wess and J. Bagger, Supersymmetry and supergravity (1992). 

32. S. Gates, M. T. Grisaru, M. Rocek, and W. Siegel, Superspace Or One Thousand and One 
Lessons in Supersymmetry, (1983), arXiv:hep-th/0108200 [hep-th]. 

33. O. A. Battistel, A. L. Mota, and M. C. Nemes, Mod. Phys. Lett. A13, 1597 (1998). 

34. O. A. Battistel, M. C. Nemes, Phys. Rev. D59, 055010 (1999), arXiv:hep-th/9811154 [hep- 
th] . 

35. D. E. Carneiro, A. Baeta Scarpelli, M. Sampaio, and M. Nemes, JHEP 0312, 044 (2003), 
arXiv:hep-th/0309188 [hep-th]. 

36. L. C. Ferreira, A. Cherchiglia, B. Hiller, M. Sampaio, and M. Nemes, Phys. Rev. D86, 
025016 (2012), arXiv:1110.6186 [hep-th]. 

37. A. Cherchiglia, M. Sampaio, and M. Nemes, Int. J. Mod. Phys. A26, 2591 (2011), 
arXiv:1008.1377 [hep-th]. 

38. E. W. Dias, A. P. Baeta Scarpelli, L. C. T. Brito, and H. G. Fargnoli, Braz. J. Phys. 40, 
228 (2010), arXiv:0912.4396 [hep-th]. 

39. C. Seijas, The Beta function of gauge theories at two loops in differential renormalization 
(2007), arXiv:0706.1793 [hep-th]. 

40. A. Vainshtein, V. I. Zakharov, and M. A. Shifman, JETP Lett. 42, 224 (1985). 

41. M. A. Shifman, A. Vainshtein, and V. 1. Zakharov, Phys. Lett. B166, 334 (1986). 

42. M. E. Peskin and D. V. Schroeder, An Introduction to quantum field theory (1995). 

43. E. Kraus, Phys. Rev. D65, 105003 (2002), arXiv:hep-ph/0110323 [hep-ph]. 

44. L. Abbott, M. T. Grisaru, and D. Zanon, Nucl. Phys. B244, 454 (1984). 

45. M. T. Grisaru, B. Milewski, and D. Zanon, Phys.Lett.B155, 357 (1985). 

46. H. Fargnoli, A. Baeta Scarpelli, L. Brito, B. Hiller, M. Sampaio, et ah, Mod. Phys. Lett. 
A26, 289 (2011), arXiv:1001.1543 [hep-th]. 



